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THREE  DIMENSIONAL  EQUILIBRIUM  AND  STABILITY  OF 
IONOSPHERIC  PLASMA  CLOUDS 


I.  INTRODUCTION 

The  dynamic  evolution  of  artificial  plasma  clouds  in  the  earth's 
ionosphere  has  been  the  subject  of  scientific  study  since  the  early  1960's 
when  the  first  barium  clouds  were  released.  The  initial  motivation  for 
these  experiments  was  of  a  diagnostic  nature;  the  motion  of  the  neutral 
and  ionized  barium  cloud  can  serve  as  a  means  to  estimate  the  magnitude  of 
ambient  electric  fields  and  neutral  winds  in  the  ionosphere.  However,  it 
was  soon  discovered  that  ionized  barium  clouds  not  only  undergo  gross 
translational  motions,  but  also  rapidly  structure  and  break  into  smaller 
pieces.  The  instability  responsible  for  this  structuring  is  the  E  x  B 

instability  which  was  first  proposed  to  explain  striations  in  weakly 

12  3 

collisional  laboratory  plasmas,  ’  and  later  applied  to  barium  clouds. 

Because  of  this  interesting  behavior,  a  considerable  amount  of 

research  has  been  devoted  to  understanding  the  dynamics  of  barium  clouds. 

The  types  of  studies  that  have  been  performed  concern  the  linear  and 

4-19 

nonlinear  evolution  of  the  E  x  B  instability,  as  well  as  large-scale 

2q  2^ 

numerical  simulations  of  ionospheric  barium  clouds.  ~  In  the  bulk  of 
these  studies,  it  has  been  assumed  that  the  ambient  magnetic  field  lines 
are  equipotentials,  i.e.,  the  electric  field  maps  perfectly  along  B,  so 
that  only  the  two  dimensional  motion  of  the  plasma  transverse  to  B  needs 
to  be  modelled.  For  the  most  part,  this  is  a  reasonable  assumption  since 
the  ratio  of  the  parallel  to  perpendicular  conductivity  in  the  F  region  of 
the  ionosphere  is  very  large  ~  10^).  However,  in  recent  studies 

this  assumption  has  been  relaxed  and  it  has  been  demonstrated  that  three 
dimensional  effects  can  have  a  dramatic  impact  on  the  gross  stability  of 
ionospheric  barium  clouds. *6-19 
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A  number 

of 

theoretical  analyses  of 

the 

influence 

of 

three 

dimensional  effects 

on 

the  evolution  of 

barium 

clouds 

have 

been 

published,16-19 

and 

3D 

numerical  simulation 

codes 

have  recently  been 

25  26  16  17 

developed.  ’  In  Sperling  et  al.  and  Drake  et  al.  the  plasma  cloud 

is  modelled  as  a  ID  slab  transverse  to  B,  and  as  a  ID  waterbag  parallel  to 

B  (the  cloud  had  a  finite  extent  along  B) .  In  this  work  the  parallel 

(along  B)  dynamics  coupled  the  unstable  modes  to  the  background  plasma 

outside  of  the  cloud.  The  distance  the  unstable  modes  map  into  the 

1/2 

background  plasma  is  roughly  given  by  X^  ~  X^  where  X^  ^ 

denotes  the  parallel  and  perpendicular  wavelength  relative  to  B, 
27 

respectively.  The  coupling  to  the  background  plasma  was  found  to  be 

stabilizing,  and  hence,  modes  with  small  transverse  wavelengths  are  more 

unstable  than  modes  with  large  transverse  wavelengths.  The  shortcomings 

of  this  model16, ^  are  that  (1)  it  is  restricted  to  modes  with  X^  «  Ln 

where  Lr  is  the  density  gradient  scale  length  of  the  plasma  cloud 

transverse  to  B  (it  should  be  noted  that  the  gross  structuring  of  plasma 

clouds  occurs  for  X^  >  Ln),  and  (2)  it  neglects  important  2D  effects 

transverse  to  B  associated  with  the  polarization  field  of  the  cloud. 

18 

These  two  shortcomings  were  redressed  in  Drake  and  Huba.  In  this 
work  a  cylindrical  2D  waterbag  cloud,  aligned  along  the  magnetic  field, 
and  immersed  in  a  uniform  neutral  wind  is  considered.  The  equilibrium 
potential  produced  by  the  neutral  wind  was  derived,  and  a  3D  stability 
analysis  was  performed.  An  important  new  physical  result  was  obtained. 
The  combination  of  finite  parallel  wavelengths  (k(|  =  2n/Xn  *  0)  and 

temperature  produces  an  azimuthal  drift  (V^)  of  the  plasma  perturbations 
of  the  cloud.  When  this  azimuthal  drift  is  greater  than  the  drift  of 
the  background  plasma  around  the  cloud  (produced  by  the  equilibrium 
potential  in  the  cloud  frame),  the  plasma  cloud  is  stable.  Physically, 


2 


this  stabilization  criterion  (V^  >  V^)  arises  because  of  the  convection  of 
modes  from  the  unstable  'backside'  of  the  cloud  to  the  stable  'frontside'. 
The  major  weakness  of  this  calculation  is  the  assumption  that  the  plasma 
cloud  is  infinite  along  B.  The  3D  effects  considered  only  arise  because 
the  perturbations  in  the  stability  analysis  have  finite  k^  The 
equilibrium  is  strictly  2D  so  that  modes  with  ^  =  0  remain  unstable. 
Although  these  results  are  significant,  direct  application  to  ionospheric 
plasma  clouds  is  tenuous  because  of  the  'arbitrariness'  of  the  parallel 
wavenumber  (k^.  The  plausible,  but  ad  hoc,  assumption  that  k(|  ~  Lz_*  ~ 
(Wrc_1’  w^ere  t^ie  half-length  of  the  cloud  along  B  and  r^  is 
its  radius  transverse  to  B,  leads  to  a  stability  criterion  which  is 

consistent  with  the  observed  size  of  barium  cloud  striations. 

19 

Subsequently,  Drake  and  Huba  completed  a  fully  3D  analysis  of 
ionospheric  plasma  clouds.  The  calculation  was  based  on  a  3D  elliptical 
waterbag  cloud  in  which  the  cloud  has  finite  extent  transverse  and 
parallel  to  B.  The  major  results  of  this  work  are  the  following.  First, 
the  equilibrium  potential  consists  of  a  polarization  potential  (caused  by 
the  neutral  wind)  and  an  ambipolar  potential  (caused  by  electron  pressure 
parallel  to  B) .  The  ambipolar  potential  causes  the  cloud  to  spin  about 
the  axis  aligned  with  the  magnetic  field.  The  rotation  rate  varies  with 
position  along  the  length  of  the  cloud  (along  B)  and  therefore  'shears' 
the  cloud.  This  shear  in  the  azimuthal  rotation  forces  the  E  x  B 
instability,  which  is  nominally  a  flute  mode  (kn  =  0),  to  have  finite  kM 
and  is  therefore  strongly  stabilizing.  Convective  stabilization  (as  was 
found  in  the  previous  cylindrical  calculation)  and  shear  stabilization  are 
additive  and  yield  a  stability  criterion  for  the  large-scale  structuring 
of  ionospheric  plasma  clouds. 
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The  results  presented  in  Ref.  19  are  restricted  to  the  special  case 
1/2 

L»z/rc  =  .  The  physical  nature  of  this  special  case  can  be 

understood  as  follows.  An  important  feature  of  our  previous  analysis  was 
the  use  of  normalized  coordinates  which  served  to  simplify  the  analysis 
and  to  demonstrate  that  the  stability  of  the  cloud  depended  only  on  a  very 
limited  number  of  parameters.  The  direction  transverse  to  B  is  normalized 
to  the  characteristic  perpendicular  scale  length  of  the  cloud  rc,  while 

*  i/2 

the  direction  parallel  to  B  is  normalized  to  L  =  r  (a., /a.)  .  In  such  a 

set  of  coordinates  the  plasma  conductivity  is  isotropic.  For  typical 

ionospheric  parameters  <r(|  >>  so  that  >>  r  .  Electrodynamically  the 

plasma  couples  much  better  along  the  magnetic  field  than  across  it  so  that 

in  this  new  coordinate  system  the  plasma  cloud  is  greatly  reduced  along  B 

in  comparison  with  the  original  physical  system.  In  these  normalized 

coordinates  the  length  of  the  cloud  along  B  is  z  *  L  /L  .  In  Ref.  19  z 

~  c  z  z  c 

=  1  so  that  the  cloud  was  a  sphere  in  our  normalized  coordinates.  In  the 
present  manuscript  we  consider  the  more  general  case  of  an  ellipsoidal 
plasma  cloud  with  z c  arbitrary.  Physically,  the  neutral  barium  clouds 
expand  radially;  subsequent  to  photo-ionization  the  plasma  barium  cloud 
continues  to  expand  along  B  but  its  motion  is  inhibited  transverse  to  B. 
Thus,  in  configuration  space,  barium  clouds  are  greatly  extended  along  B. 
They  are  ellipsoidal  in  appearance  with  the  major  axis  aligned  along  B. 
However,  in  our  normalized  coordinate  system  this  is  not  the  case.  When 
zc  <  1  the  cloud  is  an  ellipsoid  with  its  major  axis  transverse  to  B;  when 
zc  =  1  the  cloud  is  spherical;  and  when  >  1  the  cloud  is  an  ellipsoid 
with  its  major  axis  parallel  to  B.  For  typical  barium  releases  in  the  F- 
region  of  the  ionosphere  z  <  1. 
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The  organization  of  the  paper  is  as  follows.  In  the  next  section  we 

present  the  basic  equations  used  in  the  analysis.  In  Section  III  an 

analytical  equilibrium  model  is  derived.  We  derive  general  expressions 

for  the  polarization  potential  (wind  driven)  and  the  ambipolar  potential 

(electron  pressure  driven).  In  particular,  we  present  contour  plots  of 

the  ambipolar  potential,  and  associated  flows  for  z =  1/3,  1,  and  3.  A 

stability  analysis  is  then  presented  and  a  criterion  is  derived  for  the 

large-scale  structuring  of  barium  clouds.  Loosely  speaking,  we  find  that 

1/2 

clouds  with  z  <<  1  (or  L  / r  <<  (a.. /a.)  )  are  more  'stable'  than  clouds 

c  z  c  v  11  1 

1/2 

with  z  >>  1  (or  L  /r  >>  (a. ./a.)  ).  In  the  final  section  we  summarize 

c  z  c  v  II  1 

our  findings  and  apply  our  results  to  artificial  ionospheric  plasma 
clouds. 


II.  BASIC  EQUATIONS 

The  general  three-dimensional  equations  for  a  warm  plasma  cloud  in  a 

uniform  magnetic  field  B  =  B  and  a  uniform  background  neutral  wind  - 

17 

ex  (see  Fig.  1)  are  given  by 
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2 

where  <J>  is  the  electrostatic  potential,  ng  =  meve/ne  is  the  parallel 

resistivity,  v  =  v  .  +  \)  ,  D,  .  =  (v.  /S2.)cT./eB  is  the  perpendicular  ion 
J ’  e  ei  en'  li  N  in  i'  i  K  K 

diffusion  coefficient,  is  the  electron-ion  collision  frequency  and  5?^ 

and  are  the  cyclotron  and  neutral  collision  frequencies  of  species  a. 
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Equation  (1)  is  the  electron  continuity  equation  and  (1)  arises  from 

charge  neutrality  (9  •  J  =  0).  In  the  continuity  equation  the  second  term 

arises  from  the  E  x  B  drift  while  the  third  term  arises  from  the  parallel 

compression  of  the  electron  fluid.  In  the  charge  neutrality  equation  the 

first  three  terms  arise  from  the  ion  cross  field  current:  the  first  from 

the  Pedersen  conductivity,  the  second  from  diffusion  associated  with  ion- 

neutral  collisions  and  the  third  due  to  the  drag  between  the  ions  and  the 

background  neutral  wind.  The  last  term  in  (2)  again  arises  from  the 

divergence  of  the  parallel  electron  current.  In  deriving  (1)  and  (2)  we 

have  neglected  magnetic  perturbations,  taken  Tg  and  T^  to  be  spatially 

uniform,  assumed  v  /Q  ,  \>.  /Q.  <<  1  so  that  Hall  terms  can  be  discarded, 

e  e  in  i 

and  have  neglected  ion  parallel  diffusion  and  perpendicular  electron 
diffusion. 

For  typical  F-region  parameters,  the  parallel  conductivity  greatly 
exceeds  the  transverse  Pedersen  conductivity  in  the  potential  equation 
(1).  As  a  consequence  the  scale  lengths  along  the  magnetic  field  greatly 
exceed  those  in  the  transverse  direction.  It  is  therefore  convenient  to 
define  a  set  of  dimensionless  variables  in  which  the  conductivity  is 
isotropic  and  all  scale  lengths  become  comparable.  For  a  plasma  cloud  of 
perpendicular  scale  length  rc  the  appropriate  parallel  scale  length  is 

L  =  r  fa  52. /v  v.  11/Z,  (3) 

z  c  V  e  i  e  in.1  ' 

where  for  simplicity  we  take  \  and  v^n  to  be  spatially  uniform.  The 
dimensionless  equations  then  become 


9n 

at 


Vi  x  z 


Vn 


vin  a  f  94-  _  an'i  n 
S3.  3z  ln  dz  T  dzj  =  °» 


(A) 
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(5) 


where 


7  •  n7$  +  3n/3y  -  T32n/3z2  =  0, 


4>  =  c  f<#>  +  T.ln(n)/e  /Br  V 
L  l  J  c  n 


(6) 


is  the  normalized  potential,  the  dimensionless  parameter  T  is  defined  as 

T  =  c(T  +  T . )/eBV  r  ,  (7) 

e  l  n  c 

and  remaining  dimensionless  variables  are  defined  by  -»  t,  -» 

Lz  3/3z  -»  3/3z  and  n/n^  -»  n,  where  n^  is  the  density  of  the  uniform 
background. 

Because  the  conductivity  in  (5)  is  isotropic  in  our  normalized 

variables,  the  scale  lengths  along  and  across  the  magnetic  field  are  now 

comparable  (3/3z  ~  7^).  The  third  term  in  (5)  arises  from  the  current 

driven  by  the  parallel  pressure  gradient.  This  term  becomes  comparable  in 

19 

size  to  the  first  two  terms  in  (5)  when  T  >  1. 

The  parallel  compression  term  in  (4)  is  proportional  to  <<  1. 

In  our  dimensionless  equations  all  other  terms  are  of  order  unity  (for  T  - 
1)  so  that  this  term  is  small  and  can  be  neglected.  The  continuity 
equation  then  simplifies  to 

3n/3t  -  V$  x  z  •  Vn  =  0.  (8) 
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The  essential  physics  behind  the  neglect  of  the  parallel  compression  term 
in  (1)  or  (4)  is  easily  understood-  The  continuity  of  the  current  [Eq. 
(2)]  requires  that  the  divergence  of  the  parallel  current  balance  the 
divergence  of  the  Pedersen  current.  The  Pedersen  drift  is  smaller  than 
the  E  x  B  drift  [second  term  ir  (1)1  by  the  factor  v.  /Q.  .  Thus,  the 
parallel  compression  "erm  in  (1)  and  (4)  is  of  order  v^n/Q^  compared  with 
the  remaining  terms. 

Equations  (5)  and  (6)  depend  explicitly  only  on  the  single  parameter 
T.  In  addition,  there  are  two  important  "hidden"  variables:  the  ratio  of 
the  cloud  to  background  density,  M  =  n^n^;  and  the  characteristic  scale 
length  of  the  cloud  parallel  to  the  magnetic  field,  =  L^/L^. 

In  the  limit  3/3z  =  0,  (5)  and  (8)  are  independent  of  rc  so  that  the 

transverse  scale  length  is  undetermined.  However,  when  the  parallel 

dynamics  are  included,  the  parameter  T  a  rc_*  enters  the  equations.  In 

general,  we  expect  the  plasma  cloud  to  be  unstable  to  the  E  x  B 

instability  when  T  is  small.  For  a  given  cloud  configuration  (aspect- 

ratio)  we  expect  that  when  T  is  sufficiently  large  it  will  dominate  the 

remaining  terms  in  the  equations  and  either  substantially  alter  or 

stabilize  the  instability.  Thus,  for  a  given  cloud  configuration  only, 

clouds  larger  than  a  minimum  scale  size  will  be  unstable.  The  calculation 

of  this  minimum  cloud  size  was  presented  previously  for  the  simplest 
19 

configurations,  that  of  a  sphere  (an  aspect-ratio  of  one)  in  our 
dimensionless  system  of  equations.  It  was  shown  that  the  cloud  was  stable 
for  T  >  Tc(M).  In  the  present  manuscript  we  extend  this  previous  work  by 
investigating  the  equilibrium  and  stability  of  arbitrary  aspect-ratio 
ellipsoids. 
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III.  EQUILIBRIUM 


We  consider  a  waterbag  cloud  in  the  shape  of  an  ellipsoid  of  radius 
one  perpendicular  to  the  magnetic  field  and  length  z along  the  magnetic 
field  (see  Fig.  1).  The  neutral  wind  polarizes  the  cloud,  producing  a 
potential  which  causes  the  cloud  to  drift  with  a  velocity  Vq  in  the 
direction  of  the  wind.  In  addition,  an  ambipolar  potential  develops 
which  balances  the  electron  pressure  along  the  magnetic  field.  The 
ambipolar  field  causes  the  cloud  to  spin  around  its  axis  along  B. 

The  equilibrium  potentials  can  be  calculated  from  (6)  once  the 
density  profile  is  specified.  The  second  term  in  (6)  drives  the 
polarization  potential  while  the  third  drives  the  ambipolar  potential. 
Since  (6)  is  a  linear  functional  of  $  in  equilibrium,  the  total  potential 
can  be  simply  expressed  as 


*0  =  *0p  + 


0a‘ 


(9) 


We  will  evaluate  both  of  these  potentials  for  the  general  ellipsoidal 
plasma  cloud. 

It  is  convenient  to  solve  for  the  equilibrium  and  subsequently  the 
stability  of  the  cloud  in  spheroidal  coordinates,  O  h  and  (J>.  These 
coordinates  are  described  in  some  detail  in  Appendix  A.  They  are  useful 
because  surfaces  of  constant  0  where  £  is  the  "radial  variable",  are 
either  prolate  or  oblate  ellipsoids.  Thus,  in  these  coordinates  the 
equilibrium  density  n^  is  simply  given  by 

nQ(0  =  mac  -0  +  1,  (10a) 
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where 


9  i  n 

£c  =  zc/(zc  -  1)  (prolate  cloud)  (10b) 

=  zc/(l  -  zc2)1/2  (oblate  cloud)  (10c) 

is  a  function  of  the  aspect  ratio  of  the  cloud  and  H  is  the  Heaviside 
function,  and  M  is  the  ratio  of  the  cloud  density  to  background  density 
(i.e.,  nc/nb). 

The  solutions  of  the  equilibrium  equations  which  follow  will  be 
specifically  carried  out  in  prolate  coordinates.  As  discussed  in  Appendix 
A,  the  corresponding  solutions  for  oblate  coordinates  can  be  simply 
obtained  by  replacing  £,  and  everywhere  by  if.  and  i£  . 

A.  Polarization  Potential 

The  polarization  potential  is  driven  by  the  neutral  wind  and 
satisfies  the  equation 

7  •  nQV$0p  .  -  9n0/3y.  (11) 

In  the  vicinity  of  the  cloud  boundary,  the  derivatives  with  respect  to  £, 
dominate  those  with  respect  to  the  angular  variables  t|  and  ♦,  and  (11) 
simplifies  to 

—  n  ft2  -  l]  —  $  -  ^ — Z-JL-I  H  §/r  _  f  )  (12) 

9E,  n0r  )  dl  *0p  "  .  2  _  9y  V’  {1/) 

^c 

where  9f,/ay  is  given  in  (A6)  of  Appendix  A.  In  prolate  ellipsoidal 
coordinates  the  center  of  the  ellipse  corresponds  to  £  =  1  while  E,  -»  ®  at 


large  distances.  Integrating  (12)  across  the  5-function,  we  calculate  the 
jump  conditions  on  *  and  its  derivatives  across  the  boundary: 


(13a) 


3*n 

_0£  _ 
3S 


3*, 


(M  +  1) 


_ 


3£ 


=  M 


(C  -  >) 


<1  - 


h2)1'2 


sin*, 


(13b) 


where  ±  denotes  the  solution  for  5  =  £  +  £  with  £  ->  0.  Away  from  the 
Op  c 

cloud  boundary,  the  potential  satisfies  Poisson's  equation,  the  solutions 
of  which  are  given  in  Appendix  A.  Selecting  the  solutions  which  are 
evanescent  at  large  £  and  regular  as  £  -*  1,  and  matching  those  solutions 
at  £  *  £c  using  the  jump  conditions  in  (13),  we  obtain  the  polarization 
potential 


with 


♦op+<S,  ♦)  =  V0  pj(h)  sin*  q}(o/q}uc), 
*0p"(S,  1,  ♦)  =  V0  pj(n)  sin*  pJ(C)/pJ(5c), 


V0  =  M 


1 -  f(M 

-  1  L 


+  1) 


pj(tc) 


v>!' 


a 


<V]-1 


(14a) 

(14b) 


(14c) 


and  P™  and  Q™  the  associated  Legendre  polynomials  with  the  prime  denoting 
a  derivative  with  respect  to  the  argument.  The  potential  inside  the  cloud 
can  be  simplified  further  by  noting  that  (see  Appendix  A) 
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r  =  (1  -  n2)1/2(S2  -  l)1/2/Uc2  -  1)1/2 

=  -  pj(n)pj'(£;)/pj'(j;c) 


or 


V  ■  -  v-  (Ud> 

Thus,  the  cloud  moves  uniformly  with  a  velocity 

?o  -  *  x  ,f0p  -  vo  *  <15a> 

in  the  direction  of  the  neutral  wind.  In  the  limit  of  a  spherical  cloud 
(zc  =  1  and  £c  -»<»),  VQ  can  be  evaluated  explicitly, 

VQ  =  M/(M  +  3)  (15b) 

19 

as  was  found  previously.  In  the  limit  of  an  extremely  prolate  cloud,  z£ 
>>  1  and  f,c  a  1,  Vq  can  be  similarly  evaluated  as 


VQ  =  *'(M  +  D, 


(15c) 


which  is  the  velocity  of  an  infinite  cylindrical  cloud. 

Finally,  to  evaluate  the  velocity  of  an  oblate  cloud  we  make  the 
replacement  £>c  -»  i£  in  (14c)  as  follows: 


V0  =  M 


[(M  +  1) 


,1 


Pl(Uc> 


.1 


Qid^) 


l--1 


(16a) 
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For  an  extremely  oblate  cloud,  z  e  ^  «  lf  the  velocity  is  given  by 


Mz 


= 


0  "  2/lt  +  Mz 


(16b) 


Thus,  V„  -»  0  as  z  ->  0. 

0  c 

B.  Ambipolar  Potential 

The  ambipolar  potential  $Qa  is  driven  by  the  parallel  electron 
pressure  gradient  and  satisfies  the  equation 

V  •  nV*0a  =  r32n/3z2.  (17) 


Since  the  source  on  the  right  side  of  (17)  is  independent  of  the  poloidal 
angle  $,  *0a  is  similarly  independent  of  and  the  potential  *Qa  therefore 
only  produces  flows  in  the  azimuthal  direction,  i.e.,  causes  the  cloud  to 
spin  around  its  z-axis. 

In  the  vicinity  of  the  cloud  boundary,  the  derivatives  of  $Qa  with 
respect  to  E.  again  dominate  those  with  respect  to  q,  and  (17)  can  be 
approximated  as 


9_ 


Qa  = 


(18) 


where  the  subscripts  z  and  £  in  (18)  denote  a  derivative  with  respect  to 
that  variable,  and  the  functions  £zz  and  £z  have  been  evaluated  in 
Appendix  A.  Integrating  (18)  across  the  boundary  once,  we  obtain  the  jump 
condition  on  3$Qa/3?;, 


n 


3*. 


Oa 


3*. 


3£ 


-  (M  +  1) 


Oa 


3* 


Mr 


Mrs 


a_  nd  -  n2) 

3n  ^c2  -  n2  ' 


(19) 


Integrating  (18)  twice,  we  obtain  the  jump  in  the  value  of  $Qa 


*0a+  -  *0a"  "  -  F  ln(M  +  X> 


2 

11  J  ,  2 

2  z 


=  -  r  ln(M  +  l)|f 


fc*  -  ») 
k2  -  A 

(fc2  -  n2)  ’ 


(20) 


where  £  and  £  have  been  evaluated  using  the  expressions  given  in 
z  z  z 

Appendix  A. 

Away  from  the  cloud  boundary,  the  potential  again  satisfies  Poisson's 
equation.  Selecting  the  bounded  solutions  in  the  regions  £  ^  from  the 
general  solution  in  Appendix  A,  and  matching  at  £  =  £c  using  the  jump 

conditions  in  (19)  and  (20),  we  obtain  an  expression  for  the  ambipolar 
potential, 


= r  =  qtt)  >  (2ia) 

n=2  n  c 

-  pn(0 

*0^  =  r  =  anPn<*>  PTTT  ’  (21b) 

n=0  nv  c 


where 


a 


n 


anln(M  +  1) 


W<*c>l 

W  J 
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x  r(M  +  n  WiVi-i 

x  [(M  1)  Pn(^)  -  MU  J  , 


(21c) 


b  =  a  -  a  ln(M  +  1), 
n  n  n 


(21d) 
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1  _  1 

2n  +  1  r  .  2_ 

9  J  db  9  9  fi  p_(b) 

-i  t  -  tr 

c 


(21e) 


„  2n  +  i  c1  ^2  n  a  n(i  -  n2: 

Pn  =  -  — -  [  d"  *c  Pn(ri)  9b  TTT^V 

£.c  b 


( 21  f ) 


Note  that  aR  =  bn  =  0  for  n  odd  since  the  cloud  is  up-down  symmetric  about 
the  horizontal  mid-plane.  Note  also  that  in  the  expression  for  $Qg+  in 
(21a)  the  sum  over  n  extends  from  n  =  2  rather  than  n  =  0.  The  n  =  0 
contribution  is  zero  since  there  is  no  net  charge  on  the  ellipsoid  and 
thus  no  monopole  contribution  to  the  potential  as  £  -»  ®.  The  corres¬ 
ponding  solution  for  an  oblate  cloud  can  be  obtained  by  simply  replacing 
f,c  and  £  in  (21)  everywhere  by  i£c  and  i£,. 

In  general  the  potential  in  (21)  involves  an  infinite  sum  and  cannot 
be  evaluated  analytically  in  closed  form.  The  exception  is  the  limit  of  a 
spherical  cloud,  £c  -»  00 ,  for  which 


I.  ,  y.  ...  M  -  ln(M  +  1)  s 

an  -  3  ln(M  +  l)Sn0  M  +  5/2  Sn2’ 


(22a) 


.  _  M  +  (2/3)(M  +  l)ln(M  +  1)  - 

n  ’  M  +  5/2  n2 ’ 


(22b) 


in  agreement  with  previous  results. 
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The  coefficients  a  and  b  can  readily  be  evaluated  numerically  for 

n  n 

an  arbitrary  ellipsoid,  and  the  potential  can  then  be  calculated.  In 
Figs.  2a,  3a  and  4a  we  plot  the  contours  of  constant  *Qa  in  the  r-z  plane 
for  the  sphere  with  z£  =  1,  a  prolate  ellipsoid  with  =  3  and  an  oblate 
ellipsoid  with  z£  =  1/3.  The  solid  (dashed)  lines  in  Figs.  2-4  correspond 
to  positive  (negative)  values  of  the  potential.  Associated  legendre 
functions  up  to  n  =  24  have  been  retained  in  the  sums  in  (21a)  and  (21b) 
to  produce  Figs.  3  and  4.  A  surface  dipole  layer  forms  on  the  top  and 
bottom  surfaces  of  the  ellipsoids  where  the  magnetic  field  (along  z) 
intersects  the  boundary  of  the  cloud.  This  dipole  layer  produces  a 
potential  drop  at  the  boundary  which  balances  the  local  electron  parallel 
pressure  gradient.  The  dipole  layers  at  the  top  and  bottom  are  oppositely 
directed  and  therefore  asymptotically  at  large  distances  yield  a  basically 
quadrupole  potential  structure. 

In  Figs.  2b,  3b  and  4b  we  plot  the  corresponding  contours  of  the 
angular  rotation  rate 


wn  =  Vn  /r 
Oa  Oa 


(23) 


wi  th 


/dr 


=  (3*0a/3OO£/3r)  +  (3$0a/H)3h/3r,  (24) 
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the  poloidal  velocity.  The  functions  3£/3r  and  3f|/3r  are  given  in 
Appendix  A.  In  all  cases  the  angular  velocity  inside  the  cloud  is 
generally  much  smaller  than  that  outside  the  cloud.  For  a  spherical 
cloud,  the  rotation  rate  inside  the  cloud  is  uniform  so  the  cloud  rotates 
as  a  solid  body.  The  rotation  is  nonuniform  for  the  oblate  and  prolate 
spheroids.  Outside  of  the  cloud  the  fluid  rotates  positively  at  the 
midplane  and  negatively  with  higher  velocity  at  the  top  and  bottom.  For  a 
prolate  spheroid,  the  dipole  layers  are  sufficiently  far  from  the  midplane 
that  the  ambipolar  potential  in  the  midplane  is  weak  and  the  corresponding 
flows  are  small.  In  the  limit  of  extreme  elongation  the  flows  would  be 
essentially  confined  to  the  two  ends  of  the  cloud  and  the  central  portion 
of  the  cloud  would  behave  like  a  simple  2-D  infinite  cylinder. 

In  Sec.  IV  we  will  show  that  the  rotation  velocity  at  the  midplane  is 
important  in  determining  the  stability  of  3-D  plasma  clouds.  In  Figs.  2c, 
3c  and  4c  we  therefore  show  the  rotation  rate  as  a  function  of  r  at  the 
midplane.  Note  that  the  rotation  rate  in  the  cloud  is  constant  for  z£  = 
1,  falls  slightly  with  increasing  r  for  <  1  and  increases  sharply  with 
increasing  r  for  zc  >  1.  In  all  cases  the  rotation  rate  changes  sign  at 
the  boundary,  and  is  a  maximum  (in  magnitude)  just  outside  the  cloud. 

Finally,  we  emphasize  that  while  in  our  normalized  coordinates  sample 
equilibria  presented  in  Figs.  2-4  range  from  oblate  to  prolate.  In  the 
original  unnormalized  physical  coordinates,  all  of  these  clouds  would 
appear  extremely  prolate.  The  normalization  squashes  the  z  direction  by 
the  square  root  of  the  ratio  of  the  perpendicular  to  parallel 
conductivity.  The  specific  configuration  of  plasma  clouds  of  interest, 
such  as  those  resulting  from  chemical  releases  in  the  atmosphere,  will  be 
discussed  in  Sec.  V. 
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IV.  STABILITY  ANALYSIS 

In  Sec.  Ill  we  evaluated  the  equilibrium  flows  for  a  general 
ellipsoidal  plasma  cloud.  We  now  investigate  the  stability  of  such 
configurations . 

With  the  neutral  wind  blowing  from  left  to  right  in  Fig.  1,  we  would 

normally  expect  the  cloud  to  be  unstable  to  flute-like  perturbations 

located  on  the  backside  as  illustrated  in  Figs.  5a  and  5b.  The  solid 

(dashed)  lines  in  Fig.  5a  show  regions  where  the  plasma  cloud  is 

convecting  outwards  (inwards)  as  can  be  seen  in  the  midplane  slice  of  the 

cloud  in  Fig.  5b.  Along  the  direction  of  the  magnetic  field  in  Fig.  5a, 

the  phase  of  the  corrugations  on  the  surface  of  the  cloud  do  not  vary  and 

thus  the  perturbations  ire  flutelike.  In  the  case  of  a  sphere,  Drake  and 
19 

Huba  found  that  local  perturbations  propagate  in  the  azimuthal  direction 
and  that  the  rate  of  propagation  depends  strongly  on  the  polar  angle 
(latitude  above  or  below  the  midplane).  In  particular  the  azimuthal  flows 
are  small  at  the  midplane  and  increase  sharply  above  and  below  as 
illustrated  in  Fig.  5a.  As  a  consequence  the  corrugations  are  convected 
as  shown  in  Fig.  5c  and  lose  their  ilute-like  character  -  the  phase  of  the 
corrugations  varies  strongly  along  B.  The  resulting  parallel  wavelength 
has  a  dissipative  influence  on  the  instability  which  we  wish  to 
investigate. 

The  equation  for  the  plasma  density  can  be  written  in  terms  of  a 
displacement  £(h,  <fr,  t)  of  the  cloud  boundary, 

Ml,  n,  t)  =  MH(S  +  *c  -  S)  +  1.  (25) 
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Thus,  the  perturbation  of  the  density  is  given  by 

n  =  -  l  3nQ/3?;,  (.26) 

where  we  work,  in  the  frame  of  reference  of  the  moving  cloud  so  that  n^  = 
ng(£)  is  time  independent.  Since  the  plasma  flow  is  incompressible,  the 
density  perturbation  in  (26)  simply  arises  from  the  displacement  of  the 
initial  cloud  boundary  and  therefore  vanishes  away  from  the  boundary. 
Inserting  (25)  into  the  continuity  equation  (8)  and  linearizing,  we  obtain 
an  equation  for  the  time  evolution  of  t, 

ft  +  (v*  x  2  •  *s)(5*0cd+  ■  *0cn(7ri  x  z  '  +  h  x  2  •  =  °’ 

(27) 

where  $  is  the  perturbed  potential,  the  subscripts  <J>,  t|  and  £.  denote  a 
partial  derivative  with  respect  to  that  variable,  and 

*0c  =  *0  +  V  <28> 

is  the  equilibrium  potential  in  the  frame  of  the  cloud. 

We  then  linearize  the  charge  neutrality  equation  (6)  to  obtain  an 
equation  for  the  perturbed  potential  $ 

7  •  „07l  -  7  •  &,057*0c  -  <1  -  V0)(En0i;,y  .  -  0. 

(29) 
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Away  from  the  cloud  boundary,  the  second,  third  and  fourth  terms  in 

(29)  vanish  since  ng^  =  0  so  that  the  potential  satisfies  Poisson's 

equation, 

V2$  =  0.  (30) 

Thus,  the  stability  calculation  can,  in  principle,  be  completed  by  solving 

(30)  inside  and  outside  the  cloud  and  deriving  a  set  of  jump  conditions 
for  $  and  its  derivatives  across  the  boundary  to  match  the  inner  and  outer 
solutions. 

To  avoid  overly  complicating  the  equations  we  will  first  simplify  the 
problem  and  retain  only  the  dominant  contributions  to  the  equations. 
Specifically,  we  expect  that  flute-like  perturbations  with  8/ 84>  >>  3/3q 
will  be  unstable.  We  therefore  limit  the  stability  calculation  to 
perturbations  for  which  3/3<j>  >>  3/3q  >>  1.  In  this  limit  the  <|>  dependence 
of  %  and  $  can  be  represented  by  the  eikonal  $,  t,  ~  exp[i3($)J  where 

3S/3<|>  s  m(  <|>)  ,  (31) 

with  m  the  local  mode  number.  As  discussed  in  the  beginning  of  this 

section,  the  variation  of  the  equilibrium  with  q  then  forces  3/3q  t  0  and 

we  find  that,  as  in  the  case  of  the  sphere,  the  modes  are  localized  in  a 
-1/2 

region  q  -  m  around  the  midplane  of  Fig.  5  where  the  shearing  of  the 

1/2 

striations  is  weakest.  In  particular,  we  find  that  3/3q  ~  m  <<  3/3<J>  ~ 
m  and  the  mode  retains  its  basic  flute  character. 
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Away  from  the  boundary,  $  satisfies  (30),  which  in  prolate 


coordinates  (A4)  can  be  written  as 


1) 


8$  2 
3E,  "  m 


(32) 


2~  2 

where  3  $/3f|  has  been  neglected  compared  with  the  second  term  in  (32)  and 
we  have  taken  h  <<  1  for  modes  localized  around  the  midplane.  The 
solution  of  (32)  for  m  >>  1  which  are  well  behaved  as  £  -*  1  and  are 
strongly  localized  around  the  cloud  boundary  at  £  =  £c> 


$  ~  exp[±  k(£  -  Sc)/£.cl 


with  the  ±  denoting  outside  and  inside  the  cloud,  respectively,  and 

k  =  m£.c2/Uc2  -  1)  (33) 

the  local  radial  wavevector.  The  final  form  of  the  potential  perturbation 
is  given  by 


*  =  ^(h)  exp[iS(<f>)]  exp[+  k  -  €.c j / £.cJexp( vt )  ,  (34) 

where  y  is  the  local  growth  rate.  The  displacement  is  given  by  the  simple 
expression 


£(h,  4>)  =  l  exp  1  i  S(  4>)  ]  exp  (yt),  (35) 

The  evaluation  of  the  jump  conditions  on  $  at  the  cloud  boundary  is 
complicated  by  the  fact  that  the  potential  I  is  singular  at  the  cloud 
boundary.  This  can  be  seen  by  noting  that  the  term  proportional  to  f  in 
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(29)  has  three  derivatives  acting  on  n^.  Thus,  even  after  integrating 
(29)  twice  with  respect  to  we  find 

*  ~  r  hov 

which  is  singular  at  the  boundary.  This  result  should  not  be  surprising. 
The  equilibrium  potential  $qc  has  a  jump  discontinuity  at  the  boundary. 
The  displacement  of  the  equilibrium  potential  therefore  produces  a 
perturbed  potential 


I  -  5  u0c/n 

which  is  singular.  The  singularity  in  the  perturbed  potential  can  be 
eliminated  by  subtracting  this  contribution  from  $,  by  defining  a  new 
dependent  variable 


$  =  $  +  l  3$0c/3£,. 

The  continuity  equation  then  simplifies  to 

ff  +  (?♦  x  z  •  -  $0cn(?h  x  z  •  ■-  0 

and  the  charge  neutrality  equation  becomes 

V  •  nQV*  -  V  •  n0?S*oa  -  V  •  ^flQ^*Qc 

-  I1  -  ’olht),  *  rK0zz  - 0 


(36) 


(37) 


(38) 
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The  second  and  fifth  terms  in  (38)  now  cancel  to  lowest  order  (compare 
with  the  equilibrium  equation  for  $qc)  s0  that  $  remains  finite  at  the 
cloud  boundary. 

Around  the  midplane  (q  =0),  the  third  term  in  (37)  is  small  and  can 

'  2 

be  neglected  and  the  second  term  simplifies  since  V<f>  x  z  •  ^  =  ( £; 

-  l)/£  so  that 
c 


ri  *  5r*+  -  0. 


(39) 


Since  this  equation  must  be  satisfied  on  both  sides  of  the  cloud  boundary, 
*  must  be  continuous  at  the  boundary  or 

r  +  l*o  U  =  r  +  l*olv  (40) 

The  evaluation  of  the  jump  in  3$/3£;  across  the  boundary  is  more  difficult. 

2  2 

Equation  (38)  is  first  multiplied  by  ( £  -  h  )  and  integrated  across  the 

boundary.  The  details  of  the  calculation  are  presented  in  Appendix  B. 
The  result  is 


(v*  -  'Hvt)-  *  -  voK  cos+  % 

*  Mr(*v  -  Mr(Ec  -  "  °-  <41> 

2 

Note  that  the  third  term  on  the  left  side  of  (41)  is  proportional  to  h 

2  2 

and  might  be  considered  small.  It  is  multiplied  by  3  /3<J>  which  is  large 
and  it  must  therefore  be  retained.  Equation  (41)  can  now  be  reduced  to  a 
single  second  order  differential  equation  for  the  angular  dependence  of 
S(h)  by  inserting  ?  and  t  in  (34)  and  (35)  into  (41)  and  eliminating  $± 
using  (39)  and  (40), 
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32  P  ,,  ,  ,._.2_2w  .  M  +  2  **c 

— £,  +  (1  +  i8)m  ri  t,  +  1  -jTp - 7j - 

an2  Mr  *c2  -  i 


(* 


Yq  COS<t>  + 


(42) 


where 


Y0  =  m[l  -  Vq)m/(M  +  2), 


(43) 


w0(<J>)  =  m[$Jcr(h  =  0)  +  (M  +  l)*0cr(n  =  0)]/(M  +  2),  (44) 


and  5  arises  from  small  parallel  diffusion  in  the  continuity  equation. 

This  term  is  necessary  to  identify  the  bounded  eigenmodes  of  (42).  Before 

solving  the  equation,  we  can  examine  the  local  dispersion  relation  by 
2  2 

simply  ignoring  3  /3t|  in  (42), 


cos<)> 


'  1W0  + 


lm 


2  Mr 


M  +  2 


(1  +  iS)h 


(45) 


The  term  proportional  to  Yq  in  (45)  represents  the  local  growth  of 
the  E  x  B  instability  for  the  elliptical  plasma  cloud.  The  instability 
grows  on  the  backside  of  the  cloud  (cos<|>  >  0)  and  is  proportional  to  the 
relative  drift  of  the  neutral  wind  and  the  plasma  cloud.  The  second  term 
in  the  right  hand  side  of  (45)  arises  from  the  azimuthal  fluid  flow  near 
the  cloud  boundary  which  convects  the  local  perturbation  in  $  and 
therefore  produces  a  finite  frequency  w^.  These  flows  vary  with  the  angle 
f|  and  therefore  the  propagation  rate  of  local  perturbations  also  varies 
with  h*  As  discussed  earlier,  such  variation  of  the  propagation  rate 
leads  to  shearing  of  the  flute  perturbations,  which  we  expect  to  be 
stabilizing.  However,  the  angular  variation  of  the  propagation  velocity 
arising  from  the  fluid  flows  is  small  compared  with  that  arising  from  the 
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third  term  on  the  right  side  of  (45)  so  the  fluid  flows  are  simply 
evaluated  at  q  =  0. 

The  third  term  in  the  right  hand  side  of  (45)  also  causes  local 

perturbations  to  propagate.  This  term  represents  the  local  diamagnetic 

frequency  for  this  geometry.  In  an  earlier  investigation  of  the  stability 

19 

of  an  infinite  cylindrical  plasma  cloud,  Drake  and  Huba  showed  that 
perturbations  with  non-zero  parallel  wavevector  kz  have  a  diamagnetic 
frequency, 


w. 


M  r  k  2 

M  +  2  1  z  ’ 


(46) 


where  the  result  has  been  written  in  our  system  of  dimensionless  units. 
In  the  present  ellipsoidal  cloud, 


2  2  2  2  2 
3V3C  ~  l  kZ/J,  Z 
z  c 


(47) 


where  k  is  given  in  (33)  and  tz  in  (A7)  and  we  have  taken  the  limit  q  << 
1.  Substituting  (47)  into  (46)  we  obtain  the  propagation  frequency  w*, 


w, 


2  mt 
m  M  +  2 


(48) 


which  is  precisely  the  frequency  given  in  (45).  An  important  feature  of 

2 

(48)  is  that  the  frequency  varies  as  q  .  This  dependence  on  q  arises  from 

the  geometry  of  the  configuration.  The  derivative  with  respect  to  £,  is 

large,  3/3£  ~  m  but  is  orthogonal  to  z  at  the  midplane,  Vz  •  V£,  ~  £,  ~ 

2 

q.  The  q  arises  from  the  projection  of  V£  along  the  magnetic  field. 
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2 

Finally,  we  emphasize  that  is  important  even  though  it  scales  as  h 

2  -1 
because  it  is  also  proportional  to  m  .  Specifically,  a>*  ~  1  when  h  ~  m 

«  1.  This  rapid  variation  of  the  mode  propagation  velocity  with  q  causes 

shearing  of  the  striations  and  is  strongly  stabilizing,  as  we  will  now 

demonstrate. 

The  bounded  solution  to  (42)  is  given  by 

t  ~  exp[im(l  +  i S/2)  h2/2]  (49) 

and  yields  the  eigenvalue 

MF  1 

y  =  r0cos<J>  -  m  w~—  — j  ~  iwo’  (50) 

z 

c 

where  Kc  has  been  eliminated  in  favor  of  the  cloud  elongation  z£  using 

(11a).  The  growth  rate  y  in  (50)  is  valid  for  both  oblate  and  prolate 

clouds.  The  solution  for  £  in  (49)  becomes  increasingly  oscillatory  as  f| 

increases,  as  a  result  of  the  shearing  of  the  striations  depicted  in  Fig. 

5.  The  associated  damping  is  given  by  the  second  term  in  the  right  side 

of  (50).  For  a  spherical  cloud,  z£  =  1,  and  the  growth  rate  in  (50) 

19 

reduces  to  that  presented  previously.  As  the  cloud  becomes  prolate,  Z£ 
>  1  and  the  damping  rate  in  (50)  is  reduced.  The  limit  of  z£  ->  ®,  that  of 
an  infinte  cylinder,  the  damping  is  zero,  as  expected.  For  an  oblate 
cloud,  zc  <  1  and  the  damping  increases.  On  the  basis  of  these  results, 
we  would  expect  oblate  plasma  clouds  to  be  more  stable  than  prolate 
clouds. 

As  mentioned  previously  the  frequency  «q  arises  from  azimuthal  fluid 
flows  at  the  cloud  boundary  and  since  the  polarization  potential  is 
proportional  to  sin<j>  [see  (14)],  Wq  is  generally  a  strong  function  of  <p. 
In  previous  investigations  of  the  stability  of  a  cylindrical  waterbag 
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18 

cloud  to  modes  with  k  *  0,  it  was  shown  that  growing  modes  localized  an 

z 

angle  where  Wq  =  0.  At  locations  where  Wq  t  0,  perturbations  convect 
to  the  stable  front  side  of  the  plasma  cloud  where  they  are  damped.  In 
the  present  investigation  the  modes  are  not  localized  in  the  azimuthal 
angle  Nevertheless,  we  expect  that  perturbations  growing  at  locations 
where  »q  =  0  will  grow  the  finite  amplitudes  while  those  growing  to 
where  t  0  will  b.  convectively  stabilized. 

When  T  <<  1,  the  ambipolar  potential  is  small,  and  simply 

arises  from  the  flow  of  the  background  plasma  past  the  cloud.  In  this 
case  there  is  a  stagnation  point  in  the  flow  at  $  =  0  so  that  =  0  at 
that  point.  More  generally,  when  r  ~  1,  the  ambipolar  flows  which  are 
independent  of  <J>,  balance  the  polarization  flows  at  a  finite  angle  <J>. 

For  T  sufficiently  small,  <|>q  =  0  and  the  local  growth  rate  from  (50) 
is  given  by  y  -  Yq.  As  T  increases,  <t>Q  is  generally  non-zero  and  the 
local  growth  rate  where  «q(<J>q)  =  0  is  given  by 

Y  =  yq  cos*0  -  m  —  .  (51) 

z 

c 

For  T  sufficiently  large,  the  damping  in  (51)  exceeds  the  driving  term  and 
the  plasma  cloud  is  stable.  Since  Yq  “  m,  the  stability  boundary  is 
independent  of  m;  therefore,  all  azimuthal  mode  numbers  m  are  stabilized 
at  the  same  value  of  T.  The  equations  for  Tc,  the  critical  value  for  T 
for  stability,  are  given  by 

rc  =  zc2(!  -  VQ)  cos<j>0,  (52a) 

where 

«0(V  rc’  M’  zc)  =  0  (52b) 
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determines  $rt.  The  cloud  is  unstable  for  T  <  T  and  stable  for  T  >  T  . 

0  c  c 

In  Figs.  6  and  7  we  show  Tc  and  $q  versus  M  for  z  =  1,  3  and  1/3.  In  all 
cases,  ♦q  decreases  from  a  positive  to  a  negative  value  as  M  increases. 
The  decrease  of  Tc  with  M  occurs  because  Vq  approaches  1  for  large  M,  the 
cloud  velocity  approaches  the  neutral  wind  velocity.  The  driving  force  of 
the  E  x  B  instability,  which  is  proportional  to  the  slip  velocity  of  the 
neutral  wind  and  plasma  cloud,  becomes  small  [I"c  ~  1  -  Vq  in  (52a)]. 

Comparing  all  three  figures,  the  threshold  r  increases  strongly  with 
increasing  elongation  of  the  cloud.  Thus,  as  expected,  oblate  clouds  are 
more  stable  than  spherical  clouds,  which  are  more  stable  than  prolate 
clouds. 

V.  SUMMARY  AND  DISCUSSION 

The  equilibrium  and  stability  of  three  dimensional  ionospheric  plasma 

clouds  has  been  investigated.  The  motivation  for  the  present  paper  is  to 

develop  a  theoretical  model  of  plasma  clouds  for  an  arbitrary  ellipsoidal 

cross-section  (i.e.,  z£  /  1  in  Fig.  1).  This  calculation  represents  a 

19 

generalization  of  the  analysis  presented  in  Drake  and  Huba  which  is 

valid  only  for  the  spherical  case  (i.e.,  z =  1). 

In  Section  II  we  presented  the  fundamental  equations  used  in  the 
analysis  [(6)  and  (8)].  It  is  important  to  note  that  these  are  valid  in 
the  regime  3/3t  <<  <<  Q^.  These  assumptions  allow  ion  inertia  and  ion 

diffusion  parallel  to  B  to  be  neglected.  This  regime  is  typical  of  barium 
clouds  released  in  the  F  region  of  the  ionosphere.  For  example,  barium 

clouds  released  at  ~  180  km  have  3/ 3 1  ~  0.1  sec  *,  v^n  ~  1  sec-*,  and  ~ 

20  sec*. 
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In  Section  III  an  equilibrium  solution  to  (6)  and  (8)  was  derived. 

The  analysis  is  carried  out  using  a  prolate  and  oblate  ellipsoidal 

coordinate  system.  We  consider  a  waterbag  density  profile:  n  =  nc  +  n^ 

inside  the  cloud  and  n  =  n,  outside  the  cloud  where  n  ,  are  the  cloud  and 

b  c,  b 

background  densities,  respectively.  The  equilibrium  potentials  consists 
of  two  parts:  the  polarization  potential  (which  is  driven  by  the  neutral 
wind),  and  the  ambipolar  potential  (which  is  driven  by  the  parallel 
electron  pressure).  In  the  cloud  frame,  the  polarization  potential  causes 
the  background  plasma  to  flow  around  the  cloud  transverse  to  B.  The 
ambipolar  potential  produces  flows  only  in  the  azimuthal  direction,  i.e., 
it  causes  the  cloud  and  background  plasma  to  rotate  about  the  z-axis  (see 
Figs.  2-4).  The  azimuthal  rotation  rate  just  inside  and  outside  of  the 
cloud  varies  strongly  with  position  on  the  surface  of  the  cloud.  The 
variation  of  these  azimuthal  flows  causes  'shearing'  of  the  conventional 
flute  perturbations  of  the  E  x  B  instability  and  is  stabilizing. 

In  Section  IV  a  stability  analysis  of  the  ellipsoidal  waterbag 
equilibrium  was  performed.  The  calculation  is  in  the  'long  wavelength 
limit'  because  of  the  waterbag  density  profile  used,  i.e.,  X  >>  Ln  where  X 
is  the  perturbation  wavelength  and  Ln  is  the  density  gradient  scale  length 
at  the  cloud  boundary  (which  is  effectively  0).  This  limit  corresponds  to 
the  observed  structuring  scale  lengths  associated  with  barium  clouds.  The 
major  result  of  this  section  is  the  marginal  stability  criterion  embodied 
in  (52).  There  are  two  components  which  make  up  this  criterion.  First, 
(52b)  requires  <0q  =  0  where  oo^  [given  in  (44)  J  is  the  local  mode 

frequency.  For  «q  =  0  the  mode  remains  stationary  on  the  unstable 
'backside'  of  the  cloud,  and  does  not  convect  to  the  stable  'frontside'. 
Second,  (52a)  results  from  the  requirement  that  y  =  0  where  y  [given  in 
(51)]  is  the  growth  rate.  The  cloud  is  stable  for  T  =  c(Tg  +  T.,)/eBVnrc  > 
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I"  and  is  unstable  for  T  <  Tc,  where  the  solution  of  the  dispersion 
relation  in  (52)  yields  Tc  =  rc(zc,M).  In  Fig.  6,  we  plot  Tc  versus  M  for 
several  values  of  in  Fig.  7  we  plot  the  corresponding  values  of  $q. 
the  azimuthal  angle  at  which  «q  =  0  at  marginal  stability. 

Two  separate  physical  effects  combine  to  produce  the  stability 
criterion  presented  in  (52)  and  Figs.  6  and  7.  First,  the  condition  that 
the  mode  not  propagate  (<0q  =  0)  forces  the  mode  to  localize  at  a  finite 

azimuthal  angle  $q  which  is  generally  nonzero,  as  can  be  seen  in  Fig.  7. 
The  non-propagating  solution  results  when  the  diamagnetic  propagation  of 
the  mode  balances  the  propagation  caused  by  background  azimuthal  flows  and 
.  has  been  discussed  in  detail  in  Ref.  18.  The  localization  of  the  mode  at 
a  finite  value  of  <p^  weakens  the  instability  since  the  neutral  wind  is  not 
blowing  directly  across  the  density  gradient  when  4>  t  0  (see  Fig.  1). 
Complete  stabilization  of  the  mode  by  this  mechanism  requires  $  >  ji/2  so 
that  the  mode  is  no  longer  localized  on  the  backside  of  the  cloud.  This 
does  not  happen  for  the  ellipsoidal  clouds  investigated  here  since  in  Fig. 
7,  <j>Q  <  n/4. 

The  second  stabilizing  mechanism  results  from  the  "shearing"  of  the 
striations  which  is  illustrated  in  Fig.  5.  The  usual  E  x  B  instability 
causes  striations  to  form  on  the  backside  of  the  cloud  as  shown  in  Figs. 
5a  and  5b.  The  striations  are  aligned  parallel  to  the  magnetic  field  so 
that  parallel  currents,  which  are  stabilizing,  are  minimized.  In  a  3-D 
ellipsoidal  plasma  cloud  the  azimuthal  plasma  flows  associated  with  the 
ambipolar  and  polarization  potentials,  and  diamagnetic  propagation  cause 
the  striations  to  propagate  azimuthally.  However,  the  propagation  rate  is 
not  uniform  and,  in  fact,  increases  away  from  the  midplane  (z  =  0)  as 
indicated  by  the  azimuthal  arrows  in  Fig.  5a.  As  a  result,  the  striations 
begin  to  twist  or  shear  as  illustrated  in  Fig.  5c.  As  a  result,  the 
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striations  are  no  longer  flutes  but  are  forced  to  have  k t  0.  This 

effect  is  stabilizing  and  gives  rise  to  an  explicit  damping  rate  [see 

2 

(51)]  which  scales  as  T/zc  .  Generally,  both  shear  stabilization  and 
convective  stabilization  combine  to  produce  Tc  shown  in  Fig.  6.  The 
exception  is  the  set  of  points  with  ^  =  0  in  Fig.  7,  where  the 
stabilization  of  the  mode  arises  entirely  from  shear  damping. 

To  obtain  a  simple  analytic  expression  for  Tc  from  the  dispersion 
relation  in  (52)  is  difficult  because  Wq,  given  in  (44),  depends  on  the 
ambipolar  potential  which  cannot  be  calculated  analytically.  However,  for 
the  entire  range  of  M  and  z presented  in  Fig.  7,  |$y|  <  n/4  so  that  cos 
~  1  in  (52a).  The  influence  of  convection  on  the  stability  of  the  mode 
is  therefore  weak.  An  approximate  dispersion  relation  can  therefore  be 
obtained  by  taking  cos  ^  ~  1  in  (52a),  leaving 

r  /z2  =  l  -  vn.  (5: 

c  c  0  v 

The  left  side  of  (53)  represents  the  shear  damping  of  the  striations  and 
the  right  side  represents  the  driving  force  of  the  instability, 
proportional  to  the  slip  velocity  between  the  neutral  wind  (one  in  our 
normalized  units)  and  plasma  cloud  (Vq).  The  cloud  velocity  Vq,  given  in 
(14c)  and  (16a),  can  be  evaluated  explicitly  as 


=  M/(M  +  2) 

z  »  1, 
c  ’ 

(54a) 

=  M/(M  +  3) 

zc  =  1, 

(54b) 

Mz  /(2/re  +  Mz  ) 
c  c 

z  «  1. 
c 

(54c) 
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=  2/  (M  +  2). 


(55) 


For  z  »  1, 
c  ’ 


The  driving  term  reduces  to  the  usual  cylindrical  result,  as  it  should. 
For  zc  large  I"c  must  be  very  large  in  order  to  stabilize  the  cloud.  For 

zc  <<C  lf 


Tc/zc2  =  (2/ n)  /  ( 2/ ti  +  Mzc),  (56) 

and  the  Tc  is  very  small.  For  M  large  both  (55)  and  (56)  yield  ~  M-\ 
consistent  with  Fig.  6.  For  M  large,  the  relative  slip  velocity  of  the 
neutral  wind  and  plasma  cloud  is  small  so  that  weak  damping  is  sufficient 
to  stabilize  the  cloud.  Thus,  we  conclude  that  short  clouds  (zc  «  1) 
with  high  M  are  easiest  to  stabilize  while  long  clouds  (zc  »  1)  with  low 
M  are  the  most  difficult  to  stabilize. 

We  now  apply  these  results  to  the  structuring  of  F  region  barium 
releases.  It  has  been  observed  that  barium  clouds  develop  field-aligned 
striations  because  of  the  E  x  B  instability,  and  effectively  break  up  into 
smaller  'clouds'.  These  smaller  'clouds'  (or  'fingers')  also  undergo 
structuring  and  break  up  into  even  smaller  'clouds'.  This  process  is 
referred  to  as  bifurcation  and,  in  principal,  could  continue  down  to  scale 
lengths  associated  with  'dissipative'  effects  (e.g.,  electron  diffusion, 
finite  Larmor  radius  effects).  However,  observations  indicate  that  the 
striations  'freeze'  at  scale  sizes  substantially  larger  than  the 
'dissipation'  scale  lengths.  For  example,  typical  ionospheric  parameters 
for  barium  clouds  released  at  ~  180  km  are  Tg  =  T^  =0.1  eV,  B  =  0.3  G,  Vn 
=  100  m/sec,  nc  =  10°  -  10'  cm-  ,  and  n^  =  10  cm-  .  For  these  parameters 
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the  dissipation  scale  lengths  are  <  12  m  while  the  'freezing'  radius  of 

28 

clouds  is  ~  200  m.  Applying  (52)  to  ionospheric  releases,  we  compute 

the  value  of  T  based  on  the  observed  freezing  scale  radius  of  ~  200  m. 

For  the  parameters  quoted  above  we  find  that  T  =0.3.  Now,  assuming  that 

L  ~  15  km,  r  ~  200  m,  5  ~  6  x  10^  sec-^,  8.  ~  20  sec~\  \>.  ~  1  sec-\ 
z  c  e  ’  i  ’  i  ’ 

3-1  *  1/9 

v  ~  2  x  10  sec  we  find  that  L  =  r  (8  8./v  v.  )  =  250  r  =50  km. 

e  z  c  e  i  e  m  c 

Thus,  z  =  L  /L  =  1/3.  We  see  then  from  Fig.  6  that  T  ~  0.3  >  T  for  all 
c  z  z  c 

M.  A  cloud  of  radius  200  m  is  therefore  stable,  consistent  with 
observations.  Thus,  the  convective  or  shear  stabilization  of  barium 
striations  in  three  dimensions  appears  to  be  a  plausible  explanation  of 
the  freezing  phenomenon. 

Although  three  dimensional  effects  can  account  for  the  freezing  of 
striations,  it  is  not  clear  what  role  they  play  in  the  determination  of 
the  freezing  scale  size.  During  the  initial  evolution  of  an  ionospheric 
barium  cloud,  the  cloud  is  clearly  not  in  equilibrium.  In  the  presence  of 
a  neutral  wind  or  ambient  electric  field,  the  cloud  steepens  on  one  side 
and  develops  a  sharp  density  gradient.  This  'steepened'  side  can  become 
unstable  to  the  E  x  B  instability,  producing  the  observed  striations. 
During  this  initial  structuring  process,  the  waterbag  cloud  model  used  in 
this  paper  is  not  appropriate,  and  hence,  the  stability  analysis  is  not 
applicable.  This  is  not  to  say  that  3D  effects  are  not  important  during 
the  onset  of  instability,  but  rather  than  the  theory  and  results  presented 
in  this  paper  cannot  be  used.  However,  after  the  initial  structuring  of 
the  cloud,  the  smaller  clouds  or  ’fingers'  can  evolve  into  waterbag-like 
equilibrium,  and  they  are  stable  to  subsequent  structuring  because  of  3D 
effects.  In  fact,  recent  3D  numerical  simulation  results  support  this 
hypothesis . ^ 
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APPENDIX  A:  Prolate  and  Oblate  Spheroidal  Coordinates 


The  equilibrium  and  stability  of  a  prolate  (oblate)  ellipsoidal 
plasma  cloud  can  be  calculated  most  easily  in  prolate  (oblate)  spheroidal 
coordinates,  £,  q  and  4>. 

In  prolate  spheroidal  coordinates  the  "radial"  variable  1  <  £.  < 
satisfies  the  equation 


2 

z 


+ 


(Al) 


where  z,  r  and  $  are  cylindrical  coordinates.  Surfaces  of  constant  £.  are 
nested  prolate  ellipsoids  which  are  extremely  elongated  along  z  for  i,  >  1 
and  approach  spheres  as  £  ■*  °°.  The  surface  of  our  ellipsoidal  plasma 
cloud  is  defined  by  £  =  E,c,  corresponding  to  a  semi-minor  axis  of  unity 
and  a  semi-major  axis,  zm,  given  by 


z 

m 


(A2) 


Thus,  the  aspect  ratio  of  the  ellipsoidal  cloud  is  determined  by  the 
parameter  For  ^  large,  the  cloud  is  nearly  spherical,  while  for  E,Q  > 
1,  the  cloud  is  highly  elongated  along  z. 

The  angular  variable  -  1  <  q  <  1  satisfies  the  equation 


(A3) 


It  has  a  value  +1,  0,  -1  along  the  positive  z  axis,  the  midplane  (z  =  0) 
and  along  the  negative  z  axis,  respectively.  The  azimuthal  angle  <<>  is  the 
same  as  in  cylindrical  coordinates. 
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An  important  feature  of  spheroidal  coordinates  is  that  solutions  of 
Laplace's  equation  in  these  variables  are  separable, 

o2*  _  ~  *)  R  /t2  n  L,  In  .2,  3_ 

-  g.2  n2  La^  l)  as  an  (1  "  n  }  an 

*  (?TT  *  777)  $]f  ■  °- 

Thus,  the  general  solution  of  (A4)  can  be  written  as 

l,m  ’ 

x  [cm  cos(m<f»)  +  dm  sin(mij>)j,  (A5) 


where  P™  and  Q™  are  the  associated  Legendre  polynomials.  The  solution  of 
(A4)  proportional  to  Q^m(h)  has  been  eliminated  since  it  is  singular  at  t| 


-  +  1. 


The  equilibrium/stability  calculation  of  an  ellipsoidal  cloud 
requires  the  evaluation  of  a  number  of  functions  which  can  be  readily 
calculated  using  (Al)  and  (A3), 


5(?2  -  l]1/2(l  -  h2)1'2^  -  l) 
(,2  -  „2) 


sin<j>, 


42  - 1)^2  - 1)1'2 
(e2  -  12) 
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a2t  £(e2  -  i)(i  -  n2) 

3,2)K2  -  l). 

a,2  '  (52  -  n2)3 

BE/  Br  =  S^2  -  l)1/2(^c2  -  1 

)i/2(3 

-  n2)1/2/(t2  -  n2), 

3h/8r  =  -  ^  -  l)1'2^2  - 

1  -  n2)1/2/(l2  -  n2) 

an/3z  =  ^(^c2  -  i)1/2( 

i  -  < 

)/(t2-n2)1/2. 

(A8) 

(A9) 

(A10) 

(All) 


Oblate  spheroidal  coordinates,  £,  q,  $  are  similar  to  prolate 
spheroidal  coordinates.  The  analogous  Eqs.  (Al)  -  (All)  for  oblate 
spheroidal  coordinates  can  be  obtained  by  simply  replacing  £  and  £c  by  i£ 
and  i£c,  respectively,  in  these  equations.  For  example,  (Al)  transforms 
to 


2 

z 


+ 


2 

r 


S2  ♦  1 


<A12) 


for  oblate  coordinates.  In  oblate  spheroidal  coordinates  the  "radial" 
variable  £,  extends  over  the  range  0  <  £  <  00 .  Surfaces  of  constant  £, 

approach  discs  as  l  -»  0  and  spheres  as  E,  -»  ®.  The  surface  of  the  plasma 
cloud  is  again  defined  by  E  =  EQ  and  has  a  semi-major  axis  of  unity  with  a 
semi-minor  axis  (along  z)  given  by 

Zm  -  Vk2  +  1)1/2‘  (A13) 

m  c  v  c  ) 
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APPENDIX  B:  Jump  Conditions  at  the  Cloud  Boundary 

In  this  appendix  we  evaluate  the  jump  in  3$/3£  across  the  cloud 

boundary.  While  this  calculation  can  be  completed  without  approximation, 

we  avoid  presenting  excessive  algebra  by  considering  only  the  case  3/8<f>  - 
-1  1/2 

m  >>  3/3ti  ~  h  ~  m  >>  1  where  the  local  azimuthal  mode  number  m  is 

defined  in  (31).  Terms  which  are  small  in  this  ordering  are  discarded. 

2  2 

The  volume  element  in  ellipsoidal  coordinates  is  (£  -  h  )d£,dti» 

2  2 

Thus,  to  calculate  the  jump  in  3$/3£,  we  multiply  (38)  by  £  /T\  and 
integrate  across  the  boundary.  The  integral  of  the  first  term  in  (38), 
evaluated  in  ellipsoidal  coordinates,  is  given  by 


I  V-n^Vt 


C 


-  u 


2)  3 


an 


n0(-r —  +  ■  2)^t]  *• 

-  1  1  -  n  '9<r 


(Bl) 


where  I  is  shorthand  notation  for  the  integral  operator 

: .  -  ,2) 

c 


with  =  fc  ±  e  with  e  small.  Since  $  is  not  singular  at  the  boundary, 
only  the  first  term  in  (Bl)  survives, 


I  .  (5C2  -  l)2(n0»J*. 


(B2) 
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From  (33)  and  (34),  we  find  that  ^  »  m  so  that  the  expression  in  (B2)  has 
a  magnitude  m  >>  1.  Only  terms  of  this  magnitude  are  retained  as  the 
integrals  over  the  remaining  terms  in  (38)  are  carried  out.  The  integral 
of  the  second  term  in  (38)  is  given  by 


I  5  *0cE  ■ 


x  [ft  "o(«2  -  *  "o  In  l1  -  i2)  ft  ♦  +  77^2)  ?  ■rocE- 


(B3) 


In  the  equilibrium  $qc  undergoes  a  jump  at  the  cloud  boundary  so  that 
8#qc/3?;  is  singular.  The  second  and  third  terms  in  (B3)  therefore  yield 
finite  contributions.  For  3/34>  »  3/3£,  >>  1,  the  third  term  in  (B3) 
dominates.  The  angular  portion  of  *Qc?;  arises  only  ^rom  the  ambipolar 
field  which  is  independent  of  <p  so  that  (B3)  simplifies  to 


"•V&Oct 


l2  l. 

c  4 


& 


dJ;  n0*0c£. 


(B4) 


where  we  have  taken  h  <<  1.  An  expression  for  ng$oc£;  can  obtained  by 
integrating  the  equilibrium  equation  for  in  (18) 


n„$A  .  =  T  l  2/  [l,  2  -  lV 

0  0c(  z  0£  c  (  c  ) 

-  r  »2(*c2  -  1)"OE/5c2- 


(B5) 
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The  integral  in  (BA)  can  then  be  carried  out, 

I  -  F  l"'2(tc2  *  ')  V  <B6> 

-1/2  « 

Note  that  since  q  ~  m  and  3S/34>  ~  m,  the  expression  in  (B6)  has 

magnitude  m,  the  same  as  (B2).  The  integral  of  the  third  term  in  (38) 
yields 

£ 

1  V^0S  ™0c  =  fc2  "Of  d* 

x  [k  ~  ^os  is +  In  il  ~  ^os  In 

+  ^2~_  l  +  i  \  n2)  k  ^OS  9^]$0c-  (B7) 

The  first  and  third  terms  in  (B7)  are  zero,  the  first  because  n^^  =  0  on 
either  size  of  the  boundary  and  the  third  because  3$qc/3<|>  =  0  at  the 
boundary.  In  the  frame  of  the  moving  cloud,  the  polarization  potential  is 
zero  inside  the  cloud  and  at  the  boundary  while  $Qa  is  independent  of  <f>. 
The  second  term  in  (B7)  is  non-zero.  However,  since  the  second  3/3q  in 
this  term  acts  on  an  equilibrium  quantity,  and  is  therefore  of  order 
unity,  we  find 

17-&int7*n  -  m1/2  «  m. 

Oc 

This  term  is  therefore  small  compared  with  the  contributions  from  the 
first  two  terms  of  (38)  and  can  be  neglected.  The  fourth  term  in  (38) 
yields , 
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1  My  -  r«  Ms  ♦  %  lii  *  *y  y  *>o*  <BB> 

The  first  term  in  (B8)  can  be  integrated  by  parts  once  to  produce  a  term 

J  d£;  ^0£,  3£.  (V’ 

1/2 

which  is  of  order  unity.  The  second  term  in  (B8)  is  of  order  m  while 
the  third  is  of  order  m.  Thus,  only  the  last  term  need  be  retained, 

^ojy  *  "  M  cos*  V  (B9) 

where  *  *  cos<|>.  The  integral  of  the  last  term  in  (38),  can  be  written  as 

*M»-J«(\2  & 

^c 

♦  h  *  2Vz  In  KM  <B10) 

All  of  the  operators  3/3J,  in  the  parenthesis  of  (BIO)  can  be  integrated  by 

parts  so  that  they  act  on  equilibrium  functions  and  are  therefore  of  order 

1/2 

one.  Each  3/3h  in  the  parenthesis  is  of  order  m  .  The  dominant  term  in 

2  2 

(BIO)  is  therefore  the  first,  which  scales  as  3  /3t|  ~  m.  Equation  (BIO) 

simplifies  to 

:M«  ■ -"(C  -  >)W  <BI1) 

where  T\z  has  been  evaluated  using  (All).  Combining  (B2),  (B6),  (B9)  and 
(Bll),  we  obtain  the  equation  for  the  jump  in  the  slope  of  ♦  given  in 
(41). 
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Fig.  1  —  Plasma  configuration  and  geometry.  The  cloud  is  an  ellipsoidal  waterbag  (defined  by  f  =  fr)  of 
aspect  ratio  z,  (in  normalized  units).  The  cross-section  is  elliptical  in  the  xz  plane  (a)  and  circular  in  the  xy 
plane  (b),  which  is  transverse  to  B. 
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Fig.  2a  —  The  contours  of  constant  ambipolar  potential  ^  and  azimuthal  rotation  rate  <%,  =  (3  4>0 v/9r)/r  in 
(a)  and  (b),  respectively,  for  zc  =  1.0,  M  =  5.0,  and  Tc  =  1.0.  The  magnetic  field  is  upward  in  (a)  and  (b) 
and  in  corresponding  plots  in  Figs.  3  and  4.  Contours  are  equally  spaced  with  solid  (dashed)  lines  representing 
positive  (negative)  values.  The  rotation  rate  ojoa  versus  v  at  the  midplane  (z  =  0)  is  shown  in  (c).  The  plasma 
cloud  (r  <  1)  rotates  as  a  solid  body. 
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Fig.  2b  —  The  contours  of  constant  ambipolar  potential  $0,,  and  azimuthal  rotation  rate  woa  =  (0  $o 0/3r)/r  in 
(a)  and  (b),  respectively,  for  zc  =  1.0,  M  =  5.0,  and  rc  =  1.0.  The  magnetic  field  is  upward  in  (a)  and  (b) 
and  in  corresponding  plots  in  Figs.  3  and  4.  Contours  are  equally  spaced  with  solid  (dashed)  lines  representing 
positive  (negative)  values.  The  rotation  rate  w versus  v  at  the  midplane  (z  =  0)  is  shown  in  (c).  The  plasma 
cloud  (r  <  1)  rotates  as  a  solid  body. 
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Fig.  2c  —  The  contours  of  constant  ambipolar  potential  $0,,  and  azimuthal  rotation  rate  co^,  =  ( d  <l>0fl/dr)/r  in 
(a)  and  (b),  respectively,  for  zc  =  1.0,  M  =  5.0,  and  Tf  =  1.0.  The  magnetic  field  is  upward  in  (a)  and  (b) 
and  in  corresponding  plots  in  Figs.  3  and  4.  Contours  are  equally  spaced  with  solid  (dashed)  lines  representing 
positive  (negative)  values.  The  rotation  rate  versus  v  at  the  midplane  (z  =  0)  is  shown  in  (c).  The  plasma 
cloud  (r  <  1)  rotates  as  a  solid  body. 


Fig.  3a  —  Plots  similar  to  those  shown  in  Fig.  2  but  with  zc  =  3.0.  The  aximuthal  flows  are  rather  strongly 
localized  around  the  top  and  bottom  of  the  cloud. 
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Fig.  3b  —  Plots  similar  to  those  shown  in  Fig.  2  but  with  zc  =  3.0.  The  aximuthal  flows  are  rather  strongly 
localized  around  the  top  and  bottom  of  the  cloud. 


Fig.  4a  —  Plots  similar  to  those  shown  in  Figs.  2  and  3  but  with  zc  =  1/3.  The  rotation  rate  inside  the  cloud 
is  strongly  non-uniform. 


Fig.  4c  —  Plots  similar  to  those  shown  in  Figs.  2  and 
is  strongly  non-uniform. 
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Azimuthal  location  of  non-propagating  solution 
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7580 

01CY 

ATTN 

CODE 

7551 

01CY 

ATTN 

CODE 

7555 

01CY 

ATTN 

CODE 

4730/E.  MCLEAN 

01CY 

ATTN 

CODE 

4752 

01CY 

ATTN 

CODE 

4730/B .  RIPIN 

20CY 

ATTN 

CODE 

2628 

01CY 

ATTN 

CODE 

1004/P.  MANGE 

01CY 

ATTN 

CODE 

8344/M.  KAPLAN 

COMMANDER 

NAVAL  SPACE  SURVEILLANCE  SYSTEM 
DAHLGREN,  VA  22448 

01CY  ATTN  CAPT.  J.H.  BURTON 

OFFICER-IN-CHARGE 
NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK.  SILVER  SPRING,  MD  20910 
01CY  ATTN  CODE  F31 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON.  DC  20376 
01CY  ATTN  NSP-2141 

01CY  ATTN  NSSP-2722/ 

FRED  WIMBERLY 
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COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN  LABORATORY 
DAHLGREN,  VA  22448 

01CY  ATTN  CODE  DF-14/R.  BUTLER 

OFFICER  OF  NAVAL  RESEARCH 
ARLINGTON,  VA  22217 

OICY  ATTN  CODE  465 

OICY  ATTN  CODE  461 

OICY  ATTN  CODE  402 

OICY  ATTN  CODE  420 

OICY  ATTN  CODE  421 

COMMANDER 

AEROSPACE  DEFENSE  COMMAND/XPD 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB ,  CO  80912 


OICY 

ATTN 

XPDQQ 

OICY 

ATTN 

XP 

AIR  FORCE 

GEOPHYSICS  LABORATORY 

HANSCOM  AFB,  MA 

01731 

OICY 

ATTN 

OPR/HAROLD  GARDNER 

OICY 

ATTN 

LKB/ 

KENNETH  S.W.  CHAMPION 

OICY 

ATTN 

OPR/ALVA  T.  STAIR 

OICY 

ATTN 

PHD/ JURGEN  BUCHAU 

OICY 

ATTN 

PHD/ JOHN  P.  MULLEN 

AF  WEAPONS  LABORATORY 
KIRTLAND  AFB,  NM  87117 
OICY  ATTN  SUL 

OICY  ATTN  CA/ ARTHUR  H.  GUENTHER 

AFTAC 

PATRICK  AFB,  FL  32925 
OICY  ATTN  TN 


COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT- PATTERSON  AFB,  OH  45433 
OICY  ATTN  NICD/LIBRARY 

OICY  ATTN  ETDP/B.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFIN  AFB,  NY  13441 

OICY  ATTN  DOC  LIBRARY/TSLD 

OICY  ATTN  OCSE/V.  COYNE 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 
OICY  ATTN  XPFS 

SAMSO/MN 

NORTON  AFB,  CA  02409 
(MINUTEMAN) 

OICY  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  EEP/A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY,  ROOM  G-042 
WASHINGTON,  DC  20545 

OICY  ATTN  DOC  CON  FOR 

A.  LABOWITZ 

DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR 

D.  SHERWOOD 


WRIGHT  AERONAUTICAL  LABORATORIES 
WRIGHT- PATTERSON  AFB,  OH  45433-6543 
OICY  ATTN  AAAI /WADE  HUNT 
OICY  ATTN  AAAI/ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  AND  ACQ 
DEPARTMENT  CF  THE  AIR  FORCE 
WASHINGTON,  DC  20330 
OICY  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731-5000 
OICY  ATTN  J.  DEAS 

ESD/SCD-4 


EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  DOC  CON  FOR 
J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 


LIVERMORE. 

,  CA  ' 

94550 

OICY 

ATTN 

DOC  CON  FOR 
TECH  INFO  DEPT 

OICY 

nTTN 

DOC  CON  FOR 
L-389/R .  OTT 

OICY 

ATTN 

DOC  COM  FOR 
L-31/R.  HAGER 
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LOS  ALAMOS  NATIONAL  LABORATORY 
P.O.  BOX  1663 


LOS  ALAMOS,  NM 

87545 

OICY 

ATTN 

J. 

WOLCOTT 

OICY 

ATTN 

R. 

F.  TASCHEK 

OICY 

ATTN 

E. 

JONES 

OICY 

ATTN 

J. 

MALIK 

OICY 

ATTN 

R. 

JEFFRIES 

OICY 

ATTN 

J. 

ZINN 

OICY 

ATTN 

D. 

WESTERVELT 

OICY 

ATTN 

D. 

SAPPENFIELD 

LOS  ALAMOS  NATIONAL  LABORATORY 
MS  D438 

LOS  ALAMOS,  NM  87545 

OICY  ATTN  S.P.  GARY 

OICY  ATTN  J.  BOROVSKY 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  87115 


OICY 

ATTN 

W. 

BROWN 

OICY 

ATTN 

A. 

THORNBROUGH 

OICY 

ATTN 

T. 

WRIGHT 

OICY 

ATTN 

D. 

DAHLGREN 

OICY 

ATTN 

3141 

OICY 

ATTN 

SPACE  PROJ  DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  B.  MURPHEY 

OICY  ATTN  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  DC  20545 

OICY  ATTN  DR.  YO  SONG 

natl.  oceanic  &  atmospheric 

ADMINISTRATION 
ENVIRONMENTAL  RESEARCH  LABS 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

OICY  ATTN  R.  GRUBB 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 

P.O.  BOX  92957 

LOS  ANGELES,  CA  90009 


OICY 

ATTN 

I. 

GARFUNKEL 

OICY 

ATTN 

T. 

SALMI 

OICY 

ATTN 

V. 

JOSEPHSON 

OICY 

ATTN 

S. 

BOWER 

OICY 

ATTN 

D. 

OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOCIATION,  INC. 
1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 


OICY 

ATTN 

L. 

SLOAN 

OICY 

ATTN 

R. 

THOMPSON 

BERKELEY 
P.O.  BOX 
BERKELEY, 

RESEARCH 

983 

CA  94701 

ASSOCIATES 

OICY 

ATTN 

J. 

WORKMAN 

OICY 

ATTN 

C. 

PRETTIE 

OICY 

ATTN 

S. 

BRECHT 

BOEING  COMPANY,  THE 

P.O.  BOX  3707 

SEATTLE,  WA  98124 

OICY 

ATTN 

G. 

KEISTER 

OICY 

ATTN 

D. 

MURRAY 

OICY 

ATTN 

G. 

HALL 

OICY 

ATTN 

J. 

KENNEY 

CHARLES  STARK  DRAPER  LABORATORY 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  92139 

OICY  ATTN  D.B.  COX 

OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
22300  COMSAT  DRIVE 
CLARKSBURG.  MD  20871 
OICY  ATTN  G.  HYDE 


CORNELL  UNIVERSITY 

DEFT  OF  ELECTRICAL  ENGINEERING 

ITHACA .  NY  U850 

'"'ICY  ATTN  D.T.  FARLEY.  JR. 
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ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 

01CY  ATTN  H.  LOGSTON 

01CY  ATTN  SECURITY/ 

(PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  WILSHIRE  BLVD. 

SANTA  MONICA,  CA  90401 

01CY  ATTN  C.G.  GABBARD 

01CY  ATTN  R.  LELEVIER 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

01CY  ATTN  SECURITY  OFFICER 

01CY  ATTN  T.N.  DAVIS 

01CY  ATTN  NEAL  BROWN 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP- 
EASTERN  DIVISION 
77  A  STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  DICK  STEINHOF 

HSS ,  INC. 

2  ALFRED  CIRCLE 
BEDFORD.  MA  01730 

01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

01CY  ATTN  DAN  MCCLELLAND 

01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  -VA  22311 


01CY 

ATTN 

J.M. 

AEIN 

01CY 

ATTN 

ERNEST  BAUER 

01CY 

ATTN 

HANS 

WOLFARD 

01CY 

ATTN 

JOEL 

BENGSTON 

INTL  TELL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY.  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 


JAYCOR 

P.O.  BOX  85154 
11011  TORREYANA  ROAD 
SAN  DIEGO,  CA  92138 

01CY  ATTN  N.T.  GLADD 

01CY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOC  LIBRARIAN 

01CY  ATTN  THOMAS  POTEMRA 

01CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORPORATION 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET 
(P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 


01CY 

ATTN 

DASIAC 

01CY 

ATTN 

WARREN  S.  KNAPP 

01CY 

ATTN 

WILLIAM  MCNAMARA 

01CY 

ATTN 

B.  GAMBILL 

LINKABIT  CORPORATION 

10453  ROSELLE 

SAN  DIEGO.  CA  92121 

01CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 

01CY  ATTN  DEPT  60-12 

01CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
3251  HANOVER  STREET 
PALO  ALTO.  CA  94304 


01CY 

ATTN 

MARTIN  WALT/ 

DEPT  52-12 

01CY 

ATTN 

W.L.  IMHOF/ 

DEPT.  52-12 

niCY 

ATTN 

RICHARD  G.  JOHNSON 
DEPT.  52-12 

01CY 

ATTN 

J.B.  CLADIS/ 

DEFT.  52-12 
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MARTIN  MARIETTA  CORPORATION 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

01CY  ATTN  R.  HEFFNER 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  02647 


OlCY 

ATTN 

N.  HARRIS 

OlCY 

ATTN 

J.  MOULE 

OlCY 

ATTN 

GEORGE  MROZ 

OlCY 

ATTN 

V.  OLSON 

OlCY 

ATTN 

R.U.  HALPRIN 

OlCY 

ATTN 

TECHNICAL  LIBRARY 
SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA.  CA  03101 
OlCY  ATTN  P.  FISCHER 

01CY  ATTN  V.F.  CREVIER 

OlCY  ATTN  STEVEN  L.  GUTSCHE 

OlCY  ATTN  R.  BOGUSCH 

OlCY  ATTN  R.  HENDRICK 

OlCY  ATTN  RALPH  KILB 

OlCY  ATTN  DAVE  SOULE 

OlCY  ATTN  F.  FAJEN 

OlCY  ATTN  M.  SCHEIBE 

OlCY  ATTN  CONRAD  L.  LONGMIRE 

OlCY  ATTN  B.  WHITE 

OlCY  ATTN  R.  ST AG AT 

OlCY  ATTN  D.  KNEPP 

OlCY  ATTN  C.  RINO 

MISSION  RESEARCH  CORPORATION 
1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  NM  87106 

OlCY  ATTN  R.  STELLINGUERF 

OlCY  ATTN  M.  ALME 

OlCY  ATTN  L.  WRIGHT 

MITRE  CORPORATION 
UESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN.  VA  22101 

OlCY  ATTN  U.  HALL 

OlCY  ATTN  U.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD 
LOS  ANGELES.  CA  90025 

OlCY  ATTN  E.C.  FIELD.  JR 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
UNIVERSITY  PARK,  PA  16802 
01  CY  ATTN  IONOSPHERIC 
RESEARCH  LAB 


PHOTOMETRICS,  INC. 
4  ARROW  DRIVE 
WOBURN.  MA  01801 


OlCY 

ATTN 

IRVING  L.  KOFSKY 

PHYSICAL 

DYNAMICS,  INC. 

P.O.  BOX 

10367 

OAKLAND, 

CA  04610 

OlCY 

ATTN 

A.  THOMSON 

R  &  D  ASSOCIATES 

P.O.  BOX 

9695 

MARINA  DEL  REY, 

CA  90291 

OlCY 

ATTN 

FORREST  GILMORE 

OlCY 

ATTN 

V.B.  WRIGHT,  JR 

OlCY 

ATTN 

U.J.  KARZAS 

OlCY 

ATTN 

H.  ORY 

OlCY 

ATTN 

C.  MACDONALD 

OlCY 

ATTN 

BRIAN  LAMB 

OlCY 

ATTN 

MORGAN  GROVER 

RAYTHEON  CORPORATION 
528  BOSTON  POST  ROAD 
SUDBURY.  MA  01776 

OlCY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  VEST  42ND  STREET 
NEW  YORK,  NY  10036 

OlCY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
10260  CAMPUS  POINT  DRIVE 
SAN  DIEGO.  CA  92121-1522 
OlCY  ATTN  L.M.  LINSON 

OlCY  ATTN  D. A.  HAMLIN 

OlCY  ATTN  E.  FRIEMAN 

OlCY  ATTN  E .  A .  STRAi'EF 

OlCY  ATTN  C . A .  SMITH 

SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
1~1"  GnnpFIPOE  DRIVE 
mc'L£A*'.  ” '  ""'inS 

"ICY  ,-iTTN  J.  COCKAYNE 

"ICY  ATTN  E.  HYMAN 
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SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 


MENLO  PARK,  CA 

94025 

OlCY 

ATTN 

J.  CASPER 

OlCY 

ATTN 

DONALD  NEILSON 

OlCY 

ATTN 

ALAN  BURNS 

OlCY 

ATTN 

G.  SMITH 

OlCY 

ATTN 

R.  TSUNODA 

OlCY 

ATTN 

D. A.  JOHNSON 

OlCY 

ATTN 

W.G.  CHESNUT 

OlCY 

ATTN 

C.L.  RINO 

OlCY 

ATTN 

WALTER  JAYE 

OlCY 

ATTN 

J.  VICKREY 

OlCY 

ATTN 

R.L.  LEADABRAND 

OlCY 

ATTN 

G.  CARPENTER 

OlCY 

ATTN 

G.  PRICE 

OlCY 

ATTN 

R.  LIVINGSTON 

OlCY 

ATTN 

V.  GONZALES 

OlCY 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 


OlCY 

ATTN 

W.P.  BOOUIST 

TRW  DEFENSE  &  1 

SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO 

BEACH. 

CA  90278 

OlCY 

ATTN 

R.K.  PLEBUCH 

OlCY 

ATTN 

S.  ALTSCHULER 

OlCY 

ATTN 

D.  DEE 

OlCY 

ATTN 

D.  STOCKWELL, 
SNIP/ 1575 

VI SI DYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MA  01803 

01CY  ATTN  W.  REIDY 

OlCY  ATTN  J.  CARPENTER 

01CY  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH.  PA  15213 

OlCY  ATTN  M.  ZABIJSKY 


